UTTG-11-09 



Asymptotically Safe Inflation 

O ' 
(N 

Steven Weinberg* 
Theory Group, Department of Physics, University of Texas 
q ! Austin, TX, 18112 

m i 

Abstract 

q ■ Inflation is studied in the context of asymptotically safe theories of gravi- 

i-C ■ tation. Conditions are explored under which it is possible to have a long 

period of nearly exponential expansion that eventually comes to an end. 

> 

in 



On 
O 



X 



*Electronic address: weinberg@physics.utexas.edu 



1 



I. Introduction 



Decades ago it was suggested that the effective quantum field theory of 
gravitation and matter might be asymptotically safe, 1 and hence ultraviolet- 
complete. That is, the renormalization group flows might have a fixed point, 
with a finite dimensional ultraviolet critical surface of trajectories attracted 
to the fixed point at short distances. Evidence for a fixed point in the quan- 
tum theory of gravitation with or without matter has gradually accumulated 
through the use of dimensional continuation, 2 the large N approximation 3 
(where N is the number of matter fields), lattice methods, 4 the truncated 
exact renormalization group, 5 and a version of perturbation theory. 6 Re- 
cently there has also been evidence that the ultraviolet critical surface is 
finite-dimensional; it has been found that even in truncations of the exact 
renormalization group equations with more than three (and up to nine) in- 
dependent coupling parameters, the ultraviolet critical surface is just three- 
dimensional. The condition that physical parameters lie on the ultraviolet 

1 S. Weinberg, in Understanding the Fundamental Constituents of Matter, ed. A. 
Zichichi (Plenum Press, New York, 1977). 

2 S. Weinberg, in General Relativity, ed. S. W. Hawking and W. Israel (Cambridge 
University Press, 1979): 700; H. Kawai, Y. Kitazawa, & M. Ninomiya, Nucl. Phys. B 
404, 684 (1993); Nucl. Phys. B 467, 313 (1996); T. Aida & Y. Kitazawa, Nucl. Phys. B 
401, 427 (1997); M. Niedermaier, Nucl. Phys. B 673, 131 (2003) . 

3 L. Smolin, Nucl. Phys. B208, 439 (1982); R. Percacci, Phys. Rev. D 73, 041501 
(2006). 

4 J. Ambj0rn, J. Jurkewicz, & R. Loll, Phys. Rev. Lett. 93, 131301 (2004); Phys. Rev. 
Lett. 95, 171301 (2005); Phys. Rev. D72, 064014 (2005); Phys. Rev. D78, 063544 (2008); 
and in Approaches to Quantum Gravity, ed. D. Orfti (Cambridge University Press). 

5 M. Reuter, Phys. Rev. D 57, 971 (1998); M. Reuter, hep-th/9605030; D. Dou & R. 
Percacci, Class. Quant. Grav. 15, 3449 (1998); W. Souma, Prog. Theor. Phys. 102, 181 
(1999); O. Lauscher & M. Reuter, Phys. Rev. D 65, 025013 (2001); Class. Quant. Grav. 
19. 483 (2002); M. Reuter & F. Saueressig, Phys Rev. D 65, 065016 (2002); O. Lauscher 
& M. Reuter, Int. J. Mod. Phys. A 17, 993 (2002); Phys. Rev. D 66, 025026 (2002); 
M. Reuter and F. Saueressig, Phys Rev. D 66, 125001 (2002); R. Percacci & D. Perini, 
Phys. Rev. D 67, 081503 (2002); Phys. Rev. D 68, 044018 (2003); D. Perini, Nucl. Phys. 
Proc. Suppl. C 127, 185 (2004); D. F. Litim, Phys. Rev. Lett. 92, 201301 (2004); A. 
Codello & R. Percacci, Phys. Rev. Lett. 97, 221301 (2006); A. Codello, R. Percacci, & C. 
Rahmede, Int. J. Mod. Phys. A23, 143 (2008); M. Reuter & F. Saueressig, 0708.1317; P. 
F. Machado and F. Saueressig, Phys. Rev. D77, 124045 (2008); A. Codello, R. Percacci, 
& C. Rahmede, Ann. Phys. 324, 414 (2009); A. Codello & R. Percacci, 0810.0715; D. F. 
Litim 0810.3675; H. Gies & M. M. Scherer, 0901.2459; D. Benedetti, P. F. Machado, & 
F. Saueressig, 0901.2984, 0902.4630; M. Reuter & H. Weyer, 0903.2971. For a review, see 
M. Reuter and P. Saueressig, to be published [0708.1317]. 

6 M. R. Niedermaier, Phys. Rev. Lett. 103, 101303 (2009). 

7 A. Codello, R. Percacci, & C. Rahmede, Int. J. Mod. Phys. A23, 143 (2008); Ann. 
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critical surface is analogous to the condition of renormalizability in the Stan- 
dard Model, and like that condition yields a theory with a finite number of 
free parameters. 

The natural arena for applications of the idea of asymptotic safety is 
the physics of very short distances, and in particular the early universe. 8 
In Section II we show how to formulate the differential equations for the 
scale factor in a Robertson-Walker solution of the classical field equations 
in a completely general generally covariant theory of gravitation. In Section 
III we apply this result to calculate the expansion rate H for a de Sitter 
solution of the classical field equations. We are interested here in solutions 
for which H is of the same order as the scale at which the couplings are 
beginning to approach their fixed point, or larger. In this case, H turns 
out in the tree approximation to depend strongly on the ultraviolet cutoff, 
indicating a breakdown of the classical approximation. We deal with this by 
choosing an optimal cutoff, which minimizes the quantum corrections to the 
classical field equations. Section IV considers more general time-dependent 
Robertson-Walker solutions of the classical field equations with an optimal 
cutoff, and explores the circumstances under which it is possible to have an 
exponential expansion that persists for a long time but eventually comes to 
an end. An illustrative example is worked out in Section V. 

We will work with a completely general generally covariant theory of 
gravitation. (For simplicity matter will be ignored here.) The effective 

Phys. 324, 414 (2009); D. Benedetti, P. F. Machado, & F. Saueressig, 0901.2984, 0902.4630 
s The implications of asymptotic safety for cosmology have been considered by A. Bo- 
nanno and M. Reuter, Phys. Rev. D 65, 043508 (2002); Phys. Lett. B527, 9 (2002); 
M. Reuter and F. Saueressig, J. Cosm. and Astropart. Phys. 09, 012 (2005). This work 
differs from that presented here, in that they consider a severe truncation of the gravita- 
tional action, including only the cosmological constant and Einstein-Hilbert terms; they 
include matter as a perfect fluid with a constant equation of state parameter to; and they 
employ a time-dependent cutoff A. For more recent similar work that is somewhat closer 
in spirit to the present paper, see A. Bonanno and M. Reuter, J. Cosm. and Astropart. 
Phys. 0708, 024 (2007); J. Phys. Conf. Ser. 140, 012008 (2008). 
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action with an ultraviolet cutoff A takes the form 9 
I A [g] = -Jd*x v^Det^ A 4 g (A) + A 2 9l (A)R + g 2a (A)R 2 



+g 2b {A)R^R^ + A- 2 g 3a (A)R 3 + A~ 2 g 3b (A)RR^ R^ u + . . . (1) 

Here we have extracted powers of A from the conventional coupling con- 
stants, to make the coupling parameters g n (A) dimensionless. Because they 
are dimensionless, these running couplings satisfy renormalization group 
equations of the form 

A-±-g n (A) = (3 n (g(A)) . (2) 

The condition for a fixed point at g n = g ntf is that /3 n (g*) = for all n. As is 
well known, the condition for the couplings to be attracted to a fixed point 
g nif as A — > oo can be seen by considering the behavior of g n (A) when it is 
near g n *. In the case where j3 n {g) is analytic in a neighborhood of g n *, near 
this fixed point we have 

Pn{g) ->■ ^B nm (g 

m 9*m ) B nm = — , (3) 

m V og m / # 

The solution of Eq. (2) in this neighborhood is 

9n(A)^g m + J2^(^) XN (4) 

n \ 1V1 / 

9 Higher derivative theories of this sort if used in the tree approximation have long 
been known to be plagued by "ghosts"; that is, poles in propagators with residues of 
the wrong sign for unitarity. This is only if the series of operators in (1) is truncated; 
otherwise propagator denominators are not polynomials in the squared momentum, and 
there may be just one pole, or any number of poles. Even with a truncated action, because 
of the running of the couplings, there is no one Lagrangian that can be used to find the 
propagator in the tree approximation over the whole range of momenta where the various 
poles occur, and it is not ruled out that all the poles have the residues of the right sign. 
For instance, ref. 6 shows that, in a theory with only the couplings <?i, g2 a , and g2b, the 
residue of the pole in the spin 2 propagator at high mass, which had usually been supposed 
to have the wrong sign (as for instance in the work of K. S. Stelle, Phys. Rev. D 16, 953 
(1977)), in fact has a sign consistent with unitarity. More generally, Benedetti et al. in ref. 
7 point out that for any truncation or no truncation, when we l ook for a pole at a four- 
momentum p, we must take the cut-off A to be proportional to \J —p 2 , so the denominator 
of any propagator takes the form p 2 + m 2 (—p 2 ). The function m 2 (— p 2 ) is a constant at 
sufficiently low \p 2 \ , and of the form cp 2 for momenta so large that the couplings are near 
their fixed point, where c is a constant, so the equation p 2 + m 2 (~p 2 ) — for the pole 
position has no solution if — c > 1. 
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where u and A at are eigenvectors and corresponding eigenvalues of the 
matrix B nm : 

B nm u% = \ N u% . (5) 

m 

It is a physical requirement that the only eigenvectors that are allowed 
to appear in the sum in Eq. (4) are those for which the real part of the 
corresponding eigenvalues are negative, so that the couplings actually do 
approach the fixed point. The normalizations of the eigenvectors that do 
appear in Eq. (4) are free physical parameters, the only free parameters of 
the theory, except that we can adjust the over- all normalization of all the 
eigenvectors as we like by a suitable choice of the arbitrary mass scale M. 
If we choose M to make the largest of the of order unity, then M is the 
cut-off scale at which couplings are just beginning to approach their fixed 
point. 

Aside from the illustrative example considered in Section V, we will 
not carry our discussion in this paper to the point of performing numerical 
calculations, which of course would require some truncation of the series of 
terms in the action (1). Our purpose here is to lay out the general outlines of 
such a calculation, for which purpose we do not need to adopt any specific 
truncation. Our results are worked out in detail for the terms explicitly 
shown in Eq. (1), but this is only for the purposes of illustration; nothing 
in this paper assumes the neglect of higher terms. For our purposes here, 
it makes no difference whether A is regarded as a sharp ultraviolet cutoff 
on loop diagrams to be calculated using the action (1), or as a momentum 
parameter (usually called k) in a regulator term added to the action, or a 
sliding renormalization scale. 

II. Robertson Walker Solutions 

In this section we consider how to find a solution of the classical gravita- 
tional field equations for the general action (1), of the flat-space Robertson- 
Walker form 

dr 2 = dt 2 -a 2 {t)dx 2 . (6) 

It would be very complicated to derive the ten classical field equations for 
a general metric that follow from an action like (1), and then specialize 
to the case of a Robertson-Walker metric. Instead, we can much more 
easily exploit the symmetries of this metric to derive a single differential 
equation for the Hubble rate H(t) = a(t)/a(t). In showing how to derive 
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this differential equation, we will be quite general, not making any use in 
this section of the assumption of asymptotic safety. 

We can use the rotational and translational symmetries of the line ele- 
ment (6) to write the components of the variational derivatives SI A /5g^ u in 
the form 



5iA[g] 



Sgij{x) 
SIa[9] 



RW 



A 4 

— Sij a-\t) M A (t) , 



Sgio(x) 



= 



RW 



RW 



(7) 

(8) 
(9) 



the subscript RW indicating that, after taking the variational derivative, 
the metric is to be set equal to the Robertson- Walker metric defined by (6). 
(The factors A 4 /6a 2 and A 4 /2 are inserted in the definitions of M A and Af A 
for future convenience.) Also, the general covariance of the action yields the 
generalized Bianchi identity 







(10) 



By using Eqs. (7)-(9) for the Robertson-Walker metric, Eq. (10) is reduced 
to the condition: 

a 2 aM A = j t (a 3 M A ) ■ (11) 

Therefore the gravitational field equations reduce here to a single differential 
equation: 

M A (t) = , (12) 

which we see ensures the vanishing of all variational derivatives 5I A [g]/ dg^. 
This result (which holds also in the presence of spatial curvature and matter) 
is the generalization of the familiar Friedmann equation, which would apply 
if only the Einstein-Hilbert term — ^/gR/lQi^G and a vacuum energy term 
were included in the gravitational action. 

We can express M A and then A/a in terms of variational derivatives of 
the action for the Robertson- Walker metric with respect to the scale factor 
a(t). Because a(t) appears in the Robertson- Walker metric only as a factor 
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a 2 (t) in gij(x.,t), we have 



SlAlff. 



RWJ 



Sa(t) 



d 3 x 2a(t)5ij x a 3 (t) 



8gij(x.,t) 



= VA*M A (t)a 2 (t) , 



RW 



(13) 

where V is the coordinate space volume (which can be made finite by im- 
posing periodic boundary conditions.) For the flat-space Robertson- Walker 
metric (gRw)/^> the action takes the general form 



RWj 



VA' 



J dt a 3 (t)X A ( 



H(t),H(t), 



(14) 



where as usual H(t) = a(t)/a(t). Here and in Eqs. (15)— (17) below, the el- 
lipsis . . . indicates a possible dependence of Xa on second and higher deriva- 
tives of H(t). (Second and higher time derivatives do not occur in Xa if 
the integrand of the action is y 7 — Det g times an arbitrary scalar function 
of the Riemann-Christoffel curvature tensor R^ upa , including of course an 
arbitrary dependence on the curvature scalar and the Ricci tensor, but we 
do not assume that this is the case.) Comparing Eq. (13) with the result of 
a straightforward calculation of the variational derivative of the action (14) 
with respect to a(t) gives 



M A 



3X A - 3H ^ + (3H + 9H 2 ) 



dt \dH 



dH 
+ 6H 



dlA 
dH 



d fdi A \ , nTT d_ /di A \ f_ /di A \ 

dt\dHj dt 2 \dHj 



+ 



(15) 



We note that a?dM.A is a time-derivative 



a 2 aM.A 




-H + 3H< 



dT A 
dH 



+H 



dt \dHj 



(16) 



Comparing with Eq. (11), we see that M A equals the term in square brackets 
in (16), up to a possible term equal to a constant divided by a 3 (t). But the 
term in square brackets is independent of the scale of a(t), as is ftf A (t), so 
there can be no term in their difference proportional to l/a 3 (i), and thus 



M A =X A -H 



dlA 

dH 



+ (-H + 3H 2 ) dlA 



dH dt\dH ) 



(17) 
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The ten classical field equations reduce for the flat-space Robertson- Walker 
metric to the single requirement that this vanishes. 

To evaluate the terms in the action for the Robertson-Walker metric 
with no spatial curvature that are explicitly shown in Eq. (1), we note that 
for this metric R = -12H 2 - 6H and R^W™ = 36H 4 + 36H 2 H + 12H 2 . 
Using these in Eq. (1) and comparing with Eq. (14) gives 

1 A = -g (A) + K- 2 gi (K)(\2H 2 + 6H) - A- 4 g 2a (A){\2H 2 + 6H) 2 
- A- 4 g 2b {A) (36H 4 + 36H 2 H + 12# 2 ) + A~ 6 c/3 a (A)(12# 2 + 6H) 3 
+ A~V%(A)(12# 2 + 6#)(36# 4 + 36H 2 H + 12H 2 ) + ... , (18) 

where now the dots . . . denote contributions from terms not shown in (1), 
some of which involve second and higher derivatives of H. From Eq. (17), 
we then have 

M A (H, H,H,...) = -g (A) + QK- 2 gi (K)H 2 
- A~ 4 g2a(A)(2lQH 2 H - 36H 2 + 72HH^j 

- A- 4 g 2b (A)(72H 2 H - 12H 2 + 24F#) 

+ A- 6 53a (A)( - 864 H 6 + 7776H 4 H + 32A0H 2 H 2 

- 432 H 3 + 216 HH(12H 2 + 6#)) 

+ A" 6 c/ 36 (A)( - 216F 6 + 2160H 4 H + 1008 H 2 H 2 - U4H 3 

+ HH(720H 2 + A32Hfj + ... . (19) 

This is the quantity that must be set equal to zero in finding a flat-space 
Robertson-Walker solution of the classical gravitational field equations. 

III. De Sitter Solutions and Optimal Cutoff 

We can now easily find the condition for a de Sitter solution of the 
classical field equations, with 

a(t) oc e Ht , (20) 
where H is constant. Setting the quantity (19) equal to zero for H(t) = H 
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gives our condition on H: 



= N A (H)=M A (H,0,0,...) 

= - ff o(A) + 691(A) (H/Af - 864 53 a (A) (H/Af 

-216 53o (A)(#/A) 6 + ... (21) 

It is easy to find solutions of Eq. (21) that have small values of H , very 
much smaller than the scale M at which the couplings begin to approach 
their fixed points. For sufficiently small H, we can take A to be much larger 
than H, and yet small enough so that the couplings appearing as coefficients 
in (1) become independent of A, and in particular 

A 4 <7o(A) -> Pv , A 2 ffl (A) -> l/16nG N , 

where py and Gn are the conventional, A-independent, vacuum energy and 
Newton constant. Then (21) has the familiar A-independent solution 

—2 8-kGnPv 

Because of the still mysterious fact that py is observed to be much less than 
G~ 2 , this value of H is much less than G -1 / 2 , and so radiative corrections 
and higher terms in (21) can be neglected. 

We will instead be interested here in looking for solutions for which H is 
roughly of the order of the scale M at which the couplings begin to approach 
their fixed points, or larger. In this case, we face a difficult choice: How 
should we choose A? On one hand, if we choose A <C H, then we can expect 
radiative corrections to the classical result (21) to be unimportant, because 
H provides a natural infrared cutoff in loop diagrams constructed using the 
action (1). But for A <C H, the sum (21) receives increasing contributions as 
we include higher and higher terms, and whether or not the series actually 
converges, it is not useful. On the other hand, if we choose A 3> H, then 
it is reasonable to suppose that the series (21) is dominated by its lowest 
terms, but for A » H there is no reason to suppose that we can neglect 
radiative corrections to the field equations. Indeed, we can see that radiative 

10 Note that this is not the result that would be obtained by setting the derivative of 
Xa{H, 0, 0, . . .) with respect to H equal to zero. For a de Sitter metric with a(t) = exp(Ht), 
the integral over t in the action lA[g] diverges at t = 00. If we integrate only from 
t = —00 to t = 0, the integral J dta 3 (t) gives a factor 1/3H, but the derivative of 
Ia(H, 0, 0, . . -)/SH with respect to H is not zero; it equals a surface term (dXA/dH)jj, 
which again gives Eq. (21). 
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corrections to the field equations are important here, because where Eq. (21) 
is dominated by its lowest terms, it gives H a strong dependence on A. (This 
is clearest in the case where A is so large that the couplings are near their 
fixed points, in which case (21) gives H proportional to A.) The whole 
point of the renormalization group equations (2) is that physical quantities 
like H should be independent of the cutoff, but in general this is true only 
when radiative corrections are included, and since Eq. (21) gives H a strong 
dependence on A when A S> H, radiative corrections evidently can not be 
neglected. 

Ideally, we should leave A undetermined, and calculate enough of the 
radiative corrections to the field equations so that H comes out at least 
approximately independent of A. This would not be easy. Instead, we can 
try to make a judicious choice of A to minimize the radiative corrections. We 
can guess that the optimal A is roughly of the order of H, where radiative 
corrections are just beginning to be important, and the higher terms in (21) 
are just beginning to be less important. This sort of guess works quite well 
in quantum chromodynamics. The radiative corrections to a process like 
e + -e~ annihilation into jets of hadrons at an energy E are accompanied 
with powers of ]n(E/A), and to avoid large radiative corrections it is only 
necessary to take A ~ E. In this way, we can use the tree approximation to 
calculate the annihilation into, say, three jets, with the renormalization scale 
of the QCD coupling taken of order E. But in our case, radiative corrections 
are more sensitive to A, and we have to make a more careful choice of A. 

To find an optimal cutoff, we note that in principle we should find H by 
solving the full quantum corrected field equations, which give a result that 
can be schematically written as 

H tTUC = 77(A) + A77(A) , (22) 

where H(A) is defined as the solution of Eq. (21), and AH (A) represents 
the effect of radiative corrections. Instead of calculating loop graphs, we 
can get some idea of the results of such a calculation by using the tree- 
approximation field equations (21), but with A chosen at a local minimum 
of the radiative corrections to H. For such an optimal A, we have 11 

Aa77(A) = . (23) 

11 This is the weakest point in our discussion. For one thing, we do not know whether 
the condition (23) gives a local minimum or maximum of the radiative corrections. Worse, 
even if the radiative corrections are minimized, we do not know that they are small. 



10 



As already mentioned, physical quantities, including the true expansion rate 
H tTue , must be independent of A, so Eq. (23) tells us also that the expansion 
rate calculated from the classical field equations is stationary at the optimal 
cut-off 

= A^(A) . (24) 
By definition, for any A we have N\(H(A)) = 0, and by differentiating 



this with respect to A and using Eq. (24) we find that the condition for an 
optimal cutoff may be put in the form 



= a ^ n ^h) 



_ _ = A A (H(A)) + B A ((H(A)) , (25) 

H=H(A) \ / \ / 



where A A arises from the explicit dependence of N A (H) on H/A: 
A A (H) = -H^=N A {H) 

= - 12 (f) <7i(A) + 5184 0^) g 3a (A) 

+ 1296 (^j g 3b (A) + ... , (26) 

and B A comes from the running of the couplings in N A : 

B A (H) = -/3o(5(A))+6/3i( 5 (A))(F/A) 2 - 864/3 3a ( 5 (A))(^/A) 6 

-216/3 3fe ( 5 (A))(^/A) 6 + ... . (27) 

We now have two equations, (21) and (25), for the two quantities H and A, 
so it is not unreasonable to expect there to be one or more solutions, with 
both A and H roughly of order M, the only mass parameter in the theory. 



IV. Time Dependence 



The de Sitter solution found in Section II describes a universe that in- 
flates eternally. For a more realistic picture of inflation, we need a solution 
that remains close to the de Sitter solution with expansion rate near H for 
a time much longer than 1/H, but that gradually evolves away from the de 
Sitter solution, so that inflation can come to an end. (We have nothing to 
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say here about the metric before the universe enters into its de Sitter phase.) 
To find such a solution, we will consider first-order perturbations of the de 
Sitter solution, of the Robertson- Walker form (6). The expansion rate will 
take the form 

H(t) = H + 5H(t) , (28) 

with \5H(t)\ <C H. Keeping only terms in (19) of first order in 5H(t), the 
field equation A/a = becomes 



co(F, A)= + a(H, A) =2 + c 2 (H, A) =3 + . . . = , (29) 

where 

(<Wa\ 



c (H,A) = H j_ = -A A (H) , (30) 



with A\ given by Eq. (26), and 
Cl (H,A) = H 



V 8H )h 



-2lQg 2a {K){^ - 7 2(fc 6 (A) 0^ +7776 fl a a (A)^ 
+ 2160 53b (A) (^\ +... , (31) 



H 



= -72 52a (A)^ - 2452 b (A) ^ + 2592<73a(A)^j 

+ 720 536 (A) (jj +... , (32) 

and so on, with the subscript H on partial derivatives meaning that after 
taking the derivatives we set H(t) = H. Eq. (29) has an obvious solution of 
the form 

SH oc exp(£ffi) , (33) 
where £ is any root of the equation 

c (H, A) + Cl (H, A) e + c 2 (H, A) £ 2 + . . . = . (34) 
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(This is a quadratic equation in the special case in which the integrand of the 
action is ^/— Detg times an arbitrary function of the curvature tensor.) For 
positive Re £, Eq. (33) represents an instability, and the number of e-foldings 
before this instability ends the exponential expansion is ~ 1/Re £. 

We would generally expect the coefficients in Eq. (34) to be of the same 
order, in which case typical solutions for £ would be of order unity, and 
inflation would either end almost immediately (if Re £ > 0) or go on forever 
(if Re £ < 0). But there are various circumstances under which we expect 
£ to be much smaller, giving a large number of e-foldings before the end of 
inflation. 12 

1. If | Co | is much less than all the other |c„|, then Eq. (34) will have a 
solution £ ~ —cq/ci, and so much less than unity. In particular, if 
we now choose A to be the optimal cutoff described in the previous 
section, then we can use the condition (25) and Eq. (30) to write 

c (H,A) = B A (H) , (35) 

According to Eq. (27), B\(H) vanishes if the couplings are at their 
fixed point, so we can conclude that it is possible to have a long but 
not eternal period of inflation if the optimal A is large enough so that 
the couplings <? n (A) are not far from their fixed point. But there is a 
limit to how close the couplings at the optimum cutoff can be to their 
fixed point. At the fixed point, the quantities (21) and (25) are both 
functions of the single parameter H/A, and it is not likely that these 
two functions would vanish at the same value of this parameter. 

2. If the couplings are not very near their fixed point, they are sensitive 
to the free parameters of the theory that characterize the particular 
trajectory in coupling-constant space on which the couplings lie, and 
it is easy to choose these couplings to make |co| as small as we like. 
For instance, where (4) applies, all the couplings are linear in the 
normalization of the eigenvectors , the only free parameters of the 
theory. In a theory of chaotic inflation, the value of these parameters 
in any big bang containing observers may be conditioned by the re- 
quirement that Co should be small enough (and have the right sign) 
to allow the bang to become big. To be specific, in order for spatial 

12 We are concentrating here on only one mode. In all cases Eq. (34) will have more 
than one solution, and we are assuming that all modes other than the one (or several) 
with Re £ small and positive either have Re £ < or for some reason are not excited. 
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curvature not to interfere with the formation of galaxies it is neces- 
sary that the universe should expand enough during inflation so that 
whatever curvature was present at the beginning of inflation would be 
decreased enough so that the curvature term in the Friedmann equa- 
tion should not dominate over the matter term when galaxies form. 13 
As is well known, the fact that spatial curvature does not dominate 
at present requires about 60 to 70 e- foldings of inflation, 14 and the 
anthropic requirement that curvature does not interfere with galaxy 
formation is almost as restrictive. But the combination of data from 
the microwave background, baryon acoustic oscillations, and type la 
supernovae distance-redshift relations has shown 15 that (within two 
standard deviations) the fractional curvature contribution to Hq 
is in the range of —0.0178 to +0.0066. It is hard to see any anthropic 
reason for a number of e-foldings large enough to reduce the curvature 
this much. 

3. Instead of cq being anomalously small, it is possible for some or all 
of the other c n to be anomalously large, in which case again £ will be 
small and the number of e-foldings will be large. For instance, we note 
that Co unlike the other c n does not involve the couplings g 2a and g 2 b, so 
if these couplings are anomalously large, as in ref. 6, then ci, C2, etc., 
will be much larger than cq, and again we will have |£| ~ |co/ci| <C 1. 



We will now apply the above results to a classic example of higher deriva- 
tive theories of gravitation, with action limited to terms with no more than 
four spacetime derivatives: 



This theory was studied by Stelle 16 as a possible renormalizable quantum 



13 B. Freivogel, M, Kleban, M. R. Martinez, and L. Susskind, J. High Energy Phys. 
0603, 039 (2006). 

14 A. Guth, Phys. Rev. D23, 347 (1981). 

15 E. Komatsu et al, Astrophys. J. Suppl. Ser. 180, 330 (2009). 
16 K. S. Stelle, ref. 9. 



V. An Example 




+g 2b (A)R^R l 



(36) 
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theory of gravitation, and has been considered recently by Niedermaier 17 
and by Benedetti et al. 18 in connection with asymptotic safety. As is well 
known, it is possible by using the Gauss-Bonnet identity to put this action 
in the form used in refs. 17 and 18: 

I A [g] = - J d 4 xy/-Drtg~ A 4 g (A) + A 2 9l (A)R + f a (A)R 2 



+f b (A)C^R 



pupa 



where Cn Upa is the Weyl tensor, and 



, _ 92b f _ 92b 

la — 92a g ' 2 



(37) 



(38) 



For this action, Eq. (21) gives the expansion rate for a de Sitter solution 
of the field equations as 



H = A^< 7o (A)/6< ?1 (A) . 



(39) 



Instead of trying to find an optimal value of A, which minimizes radiative 
corrections to Eq. (39), here we will simply assume that A is large enough 
so that the couplings g n (A) are near their fixed point g n *, and use Eq. (39) 
to express A in terms of H: 



A = HJGgu/go* 



(40) 



with H left undetermined. 

The critical question for this sort of theory is whether the de Sitter so- 
lution has an instability that ends the eepxonential expansion after a finite 
but large number of e-foldings. As we have seen, for any small perturbation 
of the de Sitter solution, a /a is a sum of terms with the time dependence 
exp(^Ht), with £ running over the roots of Eq. (34). We are now consid- 
ering an action whose integrand is \J— Detg times a scalar function of the 
metric and the Riemann-Christoffel curvature tensor, so as remarked in the 
previous section, this equation is quadratic: 



c + ci£ + c 2 f = . 



(41) 



17 M. R. Niedermaier, ref. 6. 

1S D. Benedetti, P. F. Machado, and F. Saueressig, ref. 7. 
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For the particular action (36), the coefficients are given by 

Co = 12 g u (H/ A) 2 = 2g * (42) 

ci = 3ca = ( - 2l6g 2a * ~ 72g 2b ,) (H/A) 4 = ( - 6g 2a , - 2g 2b ,)g 2 J g 2 u , (43) 
so Eq. (41) reads 

f + 3£ = A , (44) 

where 

A = -^= , 391 T . (45) 

C 2 go* (3 52a* +92b*) 

We get a realistic picture of inflation if it turns out that A is small and 
positive. In this case Eq. (44) has a root with £ ~ —3, corresponding to a 
perturbation to a/a that decays as exp(— 3Ht), and a root with £ ~ A/3, 
corresponding to a slowly growing perturbation, that ends the exponential 
phase after about 3/ A e- foldings. 

Unfortunately, the numerical results obtained in ref. 17 and 18 are not 
encouraging. The calculations of ref. 17 are expressed in terms of coupling 
constants A, gw, u, and s, related to the couplings in Eq. (36) by 

9o = 2X/g N , gi = l/g N 

g 2a = -{l + co)/3s , g 2b = l/s. (46) 

Using a version of perturbation theory, ref. 17 found that for A — > oo the 
parameters w, A and gjy approach the fixed point values 

lu* = -0.0228 , A* = 12.69 g N */(^) 2 = 0.4227 , (47) 

while s(A) vanishes as 

s(A) -»• 11.88/ ln(A/M) , (48) 

where M is some unknown large mass. Then Eq. (45) gives 

3s 0.92 
~ ~2^g~^ ln(A/M) ' ( ' 

so A is positive, but A/M would have to be about 10 8 to give 60 e-foldings 
before inflation ends. 
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In ref. 18, by using the truncated exact renormalization group equations, 
a fixed point is found with (in our notation) 

g ^ = -0.0042 , g u = -0.0101 , g 2a * = -0.0109 , g 2b * = 0.01 . (50) 

Using these results in Eq. (45) gives A = 3.05. This is positive, but un- 
fortunately not at all small. The two roots of Eq. (44) are £ = —3.80, 
correspondign to a rapidly decaying mode, and £ = 0.80, corresponding to 
an instability that ends inflation after only a few e-foldings. 
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